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ABSTRACT 

We present the discovery of 130.87 day periodic near-infrared flux variability 
for the Class II T Tauri star WL 4 (= 2MASS J16271848-2429059, ISO-Oph 
128). Our data are from the 2MASS Calibration Point Source Working Database, 
and constitute 1580 observations in J, H and K s of a field in p Ophiuchus used 
to calibrate the 2MASS All-Sky Survey. We identify a light curve for WL 4 
with eclipse amplitudes of ~0.4 mag lasting more than one-quarter the period, 
and color variations in J-H and H-K s of ~0.1 mag. The long period cannot be 
explained by stellar rotation. We propose that WL 4 is a triple YSO system, with 
an inner binary orbital period of 130.87 days. We posulate that we are observing 
each component of the inner binary alternately being eclipsed by a circum-binary 
disk with respect to our line of sight. This system will be useful in investigating 
terrestrial zone YSO disk properties and dynamics at ~1 Myr. 

Subject headings: stars: variables: other, circumstellar matter, stars: pre-main 
sequence 



1. INTRODUCTION 



When stars form and contract onto the main sequence, remnant material can remain 
in a circumstellar disk. More than half of solar-type young stella r objects (YSOs) possess 
these primordial disks from which planets can eventually form (IMeyer et al.l 119971 ). The 
physical mechanisms responsible for the evolution and dissipation of primordial disks aren't 
directly observed, p Oph iuchus (p Oph) is a ~135 pc star- f orming region containing several 
hundred ~1 Myr YSOs (jMamajekHiooa iNatta et alibood : iBarsonv et alll2005l : lLadalll987i : 
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Lada fc Wilkingill984j ). Photometric variability is a common property of YSOs, and several 
large-sky and targeted variability studies of YSOs have been undertaken in the near-i nfrared 
(near-IR, lAlves de Oliviera II2OO8I : iBarsonv et al1ll997i : Carpenter et al.lboOlal . l20Q2h . With 
the Two Micron All- Sky Survey (2MASS) Calibration Point Source Working Database (Cal- 
PSWDB, ISkrutskie et al.ll2006l ). we are carrying out a program to study the near-IR vari- 
ability of YSOs in p Oph as a probe of stellar and circumstellar disk evolution. In this letter, 
we present the discovery of a YSO with long-term periodic variability that we attribute to 
eclipses by a circum-binary disk. 



WL 4 is a previously unremarkable ~1 Myr Clas s II T T auri star in p Oph (INatta et al. 



20061 : IStrom et all Il995f ). For WL 4, INatta et all (120061 ) estimate a J-band extinction 



of Aj=5.5, an effective st ellar temper a ture of 3715 K, an intrinsic luminosity of 1.2 L G 



0) 



and a mass of 0.45 M Q . INatta et al.l (120061 ) derive only upper limits for accretion from 



the Pa/3 and Br7 emission line equivalent widths, implying WL 4 is not a strong accre- 
tor (\og(M acc /(M Q /yr))<— 9.1). For compariso n. IStrom et al.l (11995T) d erived an extinc- 
tion Aj=5.0 and ~ 1 Myr mass of 1.4 M from iD'Antona fc Mazzitellil (119941 ) isochrones. 



Tsuboi et al 



Ratzka et al 



(2000) identify three consecutive X-ray flares with a quasi-period of ~20 hours. 
(120051 ) report a 0."176 companion from high angular resolution imaging, with 
a projected separation of ~24 AU and a flux density ratio of 0.602±0.062. Observations of 
WL 4 at infrared and sub-mm wavelengths are presented in §2.2. 



2. OBSERVATIONS 

2.1. 2MASS Calibration Observations 

2MASS imaged the entire sky in three near infrared bands between 1997 and 2001. Pho- 
tometric calibration for 2MASS was accomplished using hourly observations of 35 selected 
calibration fields, with different fields visited each hour. One of these fields in p Oph covers a 
region 8.5' wide in R.A. by 60' long in dec, centered at (R.A.,Dec) = (246.80780°,-24.68901°). 
The calibration fields were observed using the same "freeze-frame" scanning strategy used for 
the main survey that yielded a net 7.8 sec exposure on the sky per scan, with six scans of the 
field taken in alternating declination directions during each hourly calibration observation. 
In 3 ~6 month visibility windows spanning 901 days, 1582 independent scans were made of 
the field in p Oph, including 1580 detections of WL 4. The raw imaging data from each scan 
of a calibration field were reduced using the same pipeline used to proce ss the main survey . 



All source extractions from all scans were loaded into the Cal-PSWDB ( jCutri et al.ll2006l ). 
We do not identify the cause of the two missing detections for WL 4 in the Cal-PSWDB, 
but it is likely an artifact of the automated processing. 
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2.2. Thermal Infrared and Sub-Millimeter Photometry 



The Cores to Disks (c 2d) Spitzer Space T elescope Legacy program surveyed star-forming 
regions including p Oph (jEvans et al.ll2003l ). The final c2d data delivery (DR4) includes 
measured flux densities from observa tions of WL 4 with MIPS with two e pochs at 24 /im, 
one epoch each at 70 /mi and 160 /mi ( Padgett et al. 20081; Rieke et al. 2004 ). a nd two epochs 



with IRAC at 3.6, 4.5, 5.8 and 8.0 /an (lAllen et al 



m prep 



Fazio et al.ll2004h . IPadeett et al. 



(120081 ) notes that WL 4 and ISO- Oph 129 fall within an ellip tical ring of 24 / xm em ission ~1 
arc-minute in diameter centered on ISO-Oph 125 and 124. iBarsony et al.l (120051 ) observed 
WL 4 at 10.8 jum and 12.5 /im The detection at 12.5 /tm appears to be inconsistent with 
th e other measurements . Obse r vations at 70 am. wi th MIPS, and in the sub-mm as reported 



in 



Andrews fc Williams! (120071 ) ; IStanke et al.l (120061 ) , are confusion limited by emission from 



ISO-Oph 124,125, and 129. In our analysis, we treat these photometry as upper limits. 



3. ANALYSIS AND RESULTS 



3.1. Periodic Variability 

For WL 4, periodic variability is apparent in the unphased data, alternating between 
bright and faint state s. We identify a period of 130.87±0.40 days using both the Lomb 



Scar* periodogram S^iel B and the period-searching algorithm of Ipiavchan et rij 



(120081 ). The phased light curve is shown in Figure la. For JD=2450000.0, the corresponding 
phase in Figure la is 0.33. We find a peak-to-peak variation of ~0.4 mag in J, H and K s , 
with a corresponding reddening of ~0.06 mag in J-H and H-K s when WL 4 is in the bright 
state. While a ~65 day period is consistent with the 2MASS results, our best physical model 
for the system (§3.3) requires that this is an alias of the true period. The period is too long 
to be associated with a ~1 Myr YSO stellar rotation period (IRebullI 1200 ll ) . 



Super-imposed upon the long-term variability, we identify a second significant period of 
variability at 4.839±0.015 days (Figure lb). We observe this additional periodic variability in 
both the faint and bright states, with a larger amplitude in J relative to K , 9 . This variability 
is consistent with rotationally modulated stellar variability (IRebullI 120011 ). 



3.2. Spectral Energy Distribution 



The c2d IRAC observations of WL 4 were fortunate to coincide with the bright and 
faint states predicted from the Cal-PSWDB data. We model the spectral energy distribu- 
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tion (SE P) in both states. We f it by inspection the photometry with reddened PHOENIX 
NextGen ( IHauschildt et al.lll999l ) synthetic stellar spectra, and a two-temperature blackbody 
dust model. For both states of WL 4, we are able to reproduce the observations (Figure 2). 
We are able to confirm the presence of an infrared excess associated with a primordial disk. 
We summarize the model parameters in Table 1. 



3.3. Model 

The long period of the near-infrared variability necessitates a binary companion, and 
the SED indicates the presence a primordial di sk. We denot e the b inary components WL 



4a and WL4b, and the companion resolved in iRatzka et al.l (120051 ) as WL 4c. All three 
components are unresolved with 2MASS and Spitzer. The observed IRAC variability implies 
a circum-binary disk around WL 4ab, but part of the infrared excess could be associated with 
a disk around WL 4c. To explain the shape of the light curve, we postulate that a component 
of the WL 4ab binary goes into obscuration and re-emerges from behind a circum-binary 
disk every 65.44 days. The binary must be inclined with respect to the disk, and the disk 
relatively close to edge-on with respect to our line of sight. 

We can solve for the brightnesses of the three components using the total stellar lumi- 
nosity derived from the SED fit, the magnitude depth of the faint state from o ne component 



being obscured, and the flux ratio observed during a predicted bright state in IRatzka et al. 



(120051 ). We find that WL 4 is comprised of three approximately equal brightness 0.6 L 
YSOs. The symmetry between the brightnesses of WL 4a and WL 4b implies that the 
binary period is 130.87 and not 65.44 days. Our model predicts that WL 4a and WL 4b 
alternate being obscured by the circum-binary disk. The similar depth and system color 
during consecutive faint states is consistent with similar spectral types for WL 4a and WL 
4b. The color variations occur during the ingress and egress phases of the eclipses, and we 
do not model this additional complexity. The color changes may indicate a slight difference 
in spectral types for WL 4a and b, or may be due to varying scattered light flux. The esti- 
mated stellar masses imply a binary separation of 0.47 AU, or ~50 stellar radii. Finally, we 
attribute the 4.84 day periodic variability to stellar rotation modulated star-spots on WL 
4c, because this is consistent with the variability being observed in both states. 



4. DISCUSSION 



We list supporting evidence for our model in §4.1-3. 
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4.1. Derived Model Parameters and Lack of Detected Accretion Signatures 



The dust luminosities we derive are sufficient to justify the shadowing model, and the 
variability is detected from 1-8 /im. Only two parameters of our model are necessarily varied 
between the bright and faint states to reproduce the SED - the hot dust luminosity which 
is important to reproduce the IRAC channels 2-4 photometry, and the obscuration of one 
stellar component in the faint state. The change in the hot dust luminosity indicates possible 
dynamical interaction or disk "warping," or simply changing illumination/heating. The hot 
dust temperatur e is consistent wit h the stable inner orbital radius of ~2-3 times the binary 
semi-major axis (jHarringtonlll977l ) . The presence of WL 4c could account for the dynamical 
origin of a disk inclined with respect to the orbit of the inner WL 4ab binary. 

The lack of a si gnificant componen t of hot ~1000 K dust and the lack of strong ac- 
cretion signatures in iNatta et al.l (120061 ) imply the lack of a massive circum-primary disk 
around WL 4a or b. The lack of detected accretion also implies that periodically driven 
accretion by the companion is no t a favored scenario , such as is hypothesized for DQ Tau 
and AA Tau (IMathieu et al.lll997l ; iBouvier et al.ll2003l ). Finally, the apparent K s magnitude 
of WL 4 is consistent with a directly vis ible p Oph YSO, and is too bright to be an edge-on 
disk system seen only in scattered light (IStapelfeldt et al.ll2008l ; IWatson fc Stapelfeldtl 120071 ; 
Stapelfeldt et al.lll997h . 



4.2. Duration of Faint State 



A stable circum-st ellar disk around WL 4a or WL 4b with an outer radius of one-third 
the binary separation (lArtymowicz fc Lubowill994j ) would produce an eclipse only ~10% of 
the period, so this configuration can be ruled out. Excepting for the ~4.84 day periodic 
variability, the light curves in the faint and bright states are relatively flat and smoothly 
varying for ~25 and ~13 days respectively. The transition between bright and faint states, 
including the "kinks" in the light curve at phases of 0.14, 0.36, 0.63 and 0.86, last ~13 days 
apiece. We speculate that the "kinks" could be due to disk substructure such as a puffed up 
or over-dense inner edge. 



4.3. Similarities to KH-15D 



KH-15D is a unique K6-7 binary YSO in NGC 2264 flKearns fc Herbstlll998l : iBadalian & Erastova 
19701 ). Every ~48 days, KH-15D periodically varies in brightness by ~4 magnitudes be- 
tween relatively quiescent bright and faint states. An inclined binary with a ~48 day pe- 
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riod surrounded by a nearly edge-on circum-binary disk is invoked to explain these obser- 



disk ( 


Winn et al. 


2006 


Johnson & Winn 


2004 



Johnson et al.ll2005t IWinn et al.ll2004l ; IChiang fc Murray- Clayl 12004 



Winn et al. 



20031 ). The near-IR photometry and colors vary in tan- 



dem in a similar fashion for both KH-15D and WL 4, although the amplitu de differs. When 



both s ources get brighter, bo th sources also get redder in J-H and H-K s (IKusakabe et al. 



20051 ). IKusakabe et al.l (120051 ) attribute the reddening for KH-15D to changes in the scat- 



tered light flux. Similar features such as the "kinks" in the light curve in Figure la have also 
evolved with time for KH-15D. 



5. CONCLUSIONS AND FUTURE WORK 

We have identified periodic variability for the YSO WL 4 in p Oph that is likely due 
alternating eclipses of two different components of a binary system by a circum-binary disk. 
The faint state lasts longer than the bright state, ruling out a co-planar stable circum-primary 
disk to explain the observations. WL 4 is a unique and valuable probe of YSO terrestrial 
zone disk evolution. We want to confirm our binary model for the light curve modulation. 
The radial velocity amplitude implied by our model should be measurable through high- 
resolution near-IR spectroscopy with adaptive optics. Detailed modeling of the disk with 
orbital dynamics of the binary is warranted to reproduce the observed light curves, including 
the "kinks" , and to investigate the dynamical stability of our model. Near- and mid-infrared 
photometric and spectroscopic observations will enable a characterization of the dust grain 
properties and disk structure. 

This publication makes extensive use of data products from 2MASS, which is a joint 
project of the University of Massachusetts and IPAC/Caltech, funded by NASA and NSF. 
This research has made use of the NASA/ IPAC Infrared Science Archive, which is operated 
by JPL, Caltech, under contract with NASA. Thanks to Mike Meyer, Mike Werner, Angelle 
Tanner, and Eric Agol for their conversations and comments. Parts of the research described 
in this publication was carried out at JPL. 



REFERENCES 

Allen, L., et al., ApJ, in prep. 

Alves de Oliviera, O, & Casali, M., 2008, A&A, in press 



- 7- 

Andrews, S., & Williams, J., 2007, ApJ, 671, 1800 

Artymowicz, P., & Lubow, S., 1994, ApJ, 421, 651 

Badalian, H., & Erastova, L., 1970, Astronomicheskij Tsirkulyar, 591, 4 

Barsony, M., Ressler, M., & Marsh, K., 2005, ApJ, 630, 381 

Barsony, M., Kenyon, S., Lada, E., & Teuben, P., 1997, ApJS, 112, 109 

Becklin, E., Matthews, K., Neugebauer, G., & Willner, S., 1978, ApJ, 220, 831 

Bouvier et al, 2003, A&A, 409, 169 

Burrows, A., et al, 2001, Rev. Mod. Phys., 73, 719 

Carpenter, J., Hillenbrand, L., Skrutskie, M., & Meyer, M., 2002, AJ, 124, 1001 
Carpenter, J., Hillenbrand, L., & Skrutskie, M., 2001, AJ, 121, 3160 
Chiang, E. I., & Murray-Clay, R. A. 2004, ApJ, 607, 913 
Chiang, E. I., et al., 2001, ApJ, 547, 1077 

Cutri, R., et al., 2006, Explanatory Supplement to the 2MASS All Sky Data Release and Ex- 
tended Mission Products, |http : / / www. ipac .caltech.edu/ 2mass / releases / allsky /do c/ 

D'Antona, F., & Mazzitelli, I., 1994, ApJS, 90, 467 

Evans, N.J., et al., 2003, PASP, 115, 965 

Fazio, G. C, et al. 2004, ApJS, 154, 39 

Harrington, R.S., 1977, A J, 82, 753 

Hauschildt, P., Allard, F., & Baron, E. 1999, ApJ, 512, 377 

Hut, P., 1981, A&A, 99, 126 

Johnson, J. A., et al. 2005, AJ, 129, 1978 

Johnson, J. A., & Winn, J. N. 2004, A J, 127, 2344 

Jura, M., Malkan, M., White, R., Telesco, C, Pina, R., & Fisher, R. S. 1998, ApJ, 505, 897 
Kearns, K., & Herbst, W., 1998, AJ, 116, 261 



- 8- 



Kusakabe, N., et al., 2005, ApJ, 632, L139 

Lada, C. J. 1987, in IAU Symp. 115, Star Forming Regions, ed. M. Peimbert & J. Jugaku ( 
Dordrecht: Reidel ), 1 

Lada, C. J., & Wilking, B. A. 1984, ApJ, 287, 610 

Mamajek, E., 2008, AN, 329, 10 

Mathieu et al, 1997, A J, 113, 184 

Mathis, J., 1990, AR&A, 28,37 

Meyer et al, 1997, A J, 114, 288 

Natta, A., Testi, L., & Randich, S., 2006, A&A, 452, 245 

Padgett, D., et al, 2008, ApJ, 672, 1013 

Plavchan, P. et al, 2008, ApJS, 175, 191 

Ratzka, T., Khler, R., & Leinert, Ch., 2005, A&A, 437, 611 

Rebull, L., 2001, AJ, 121, 1676 

Rieke, G., et al, 2004, ApJS, 154,25 

Scargle, J., 1982, ApJ, 263, 835 

Siess, L., Dufour, E. & Forestini, M., 2000, A&A, 358, 593 
Skrutskie, M., et al., 2006, AJ, 131, 1163 
Stanke, T., et al., 2006, A&A, 447, 609 
Stapelfeldt, K., et al., 2008, in prep. 

Stapelfeldt, K., et al. 1997, in Proceeding of IAU Symposium 182, B. Reipurth and C. 
Bertout editors, Kluwer Academic Publishers, pp 355-364 

Strom, K., Kepner, J., & Strom, S., 1995, ApJ, 438, 813 

Tsuboi, Y., et al., 2000, Adv. Space Res., 25, 314, 535 

Watson, A. & Stapelfeldt, K., 2007, A J, 133, 845 

Winn, J.N., et al., 2006,ApJ, 644, 510 



Winn, J.N., Holman, M.J., Johnson, J.A., Stanek, K.Z., & Garnavich, P.M., 2004, ApJ, 603, 
L45 

Winn, J. N., Garnavich, P. M., Stanek, K. Z., & Sasselov, D. D. 2003, ApJ, 593, L121 



This preprint was prepared with the AAS IATgX macros v5.2. 



-10- 



Table 1. SED Model Parameters 



Parameter 




Value 


Fixed 


Distance 


135pc 




Short A extinction power law 


-1.7 a 




Long A extinction power law 


-l a 




Extinction power law transition 


3.5 /im' 


i 


T D ,T C , R 6 ,R C , M b , M c 


T a , Ra 


, Ma b 


Varying - Best Fit 


WL 4 Aj extinction 


6 mag 




T a , Ra, M a 


3600 K 


c , 2.0 Rq, 0.40 M Q d 


cold dust T, L 


160 K e 


, 0.017 L 


hot dust T, L bright state 


560 K e 


, 0.055 L 


hot dust T, L faint state 


560 K e 


, 0.061 L 



"Ex t inction wavelength dependence adopted from iMathisI 
il990h ; lBecklin et alj <197Sh . 

b An equal-mass triple system model is used (§3.3). 

c Fit to nearest 100 K. 

d E stimated from 1 Myr 3600K Seiss isochrone. llSeiss et al.l 

booch . 

e Fit to nearest 5 K. 
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Fig. 1.— (a) Left: Top Panel - K s band Cal-PSWDB light curve data in black for WL 4. 
Middle Panel: J-H Cal-PSWDB color curve. Bottom Panel: H-K s Cal-PSWDB color curve. 
Data are folded to a period of 130.87 days, and plotted as a function of period phase, (b) 
Right: Portions of the J-band light curve data during the faint state to show the star-spot 
variability modulated by the ~5 day stellar rotation period. Overlaid in red are the best 
fitting sinusoids plus linear ramps, which reduce the photometric scatter of ~0.07 mag by 
a factor of ~2. The residuals are structured and indicate the linear ramps are an over- 
simplification; we do not identify the physical origin of these variations. The sinusoids have 
different phases for all of the spans of faint or bright states that we investigate, implying star- 
spot evolution with a coherence time-scale longer than the rotation period and on the order 
of the binary period. Sinusoidal variations are largest in amplitude at J-band relative to 
K s -band. For all panels, each group of six scans from a single hourly calibration observation 
are co-added, and 1-a error bars are shown in teal. 



-12- 




Fig. 2.— Spectral Energy Distribution for WL 4. DATA: Red circles - 2MASS Cal-PSWDB 
flux densities in the faint state (Figure la phase 0); blue circles - 2MASS Cal-PSWDB flux 
densities in the bright state (phase 0.25); red squares - c2d IRAC flux densities in the faint 
state (JD=2453071.831, Figure la phase 0.142); blue squares - c2d IRAC flux densities in 
the bright state (2453092.736, 0.302); green diamonds - 10 .8 and 12.5/im measurements 
(-2450198.0, 0.183; -2451358.0, 0.047. iBarsony et all 120051 ); purple diamond - c2d MIPS 
24 /an flux density at both epochs (2453083.888, 0.235; 2453084.113, 0.236); purple triangle 

- c2d MIPS 70 jum flux density upper limit; red triangles - sub-mm flux density upper limits 
(lAndrews fc Williamsll2~0~0~7l ; Istanke et alfcood ) . The IRAC epochs correspond to the channel 
1 and 3 observations. The channel 2 and 4 observations were taken within one hour of these 
times. MODEL: solid red line - cold dust component; solid orange line - hot dust component 
in the faint state; dashed orange line - hot dust component in the bright state; cyan line 

- reddened synthetic spectra for the two WL 4 components un-obscured in the faint state; 
magenta line - reddened synthetic spectra for the three WL 4 components un-obscured in 
the bright state; black solid line - sum of stellar components, hot dust and cold dust spectra 
in faint state; black dashed line - sum of stellar components, hot dust and cold dust spectra 
in bright state. 



